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ABSTRACT 

This  paper  presents  two  distributed  algorithms  for  finding 
shortest  paths  from  a  source  node  to  all  other  nodes  in  an  N-node 
network.  These  algorithms  are  executed  at  individual  nodes  using 
only  local  information.  Algorithm  1  works  in  networks  where  there 
are  no  topological  changes  such  as  link  failures,  link  recoveries 
or  changes  of  link  lengths.  Algorithm  2  is  a  modification  of 
Algorithm  1  for  networks  where  there  are  topological  changes. 

Algorithm  1  determines  the  optimal  shortest  paths  in  at  most 

3 
N  /4   steps,  which  is  only  one-half  of  the  computational  upper 

bounds  of  Abram  and  Rhodes'  and  Segall,  Merlin  and  Gallager's 

algorithms.   After  the  last  topological  change,  Algorithm  2 

determines  the  optimal  shortest  paths  in  the  same  number  of  steps 

as  Algorithm  1.   There  are  many  situations  where  the  present 

algorithms  will  work  up  to   N/2   times  faster  than  the  algorithms 

proposed  by  these  authors. 


1.   INTRODUCTION 

Finding  shortest  paths  in  a  network  is  an  important 
problem  and  has  many  applications  in  computer  networks  [2], 

[5],  [10].   Two  types  of  algorithms  are  available  for  finding 
the  shortest  paths  in  networks.   The  first  type  of  algorithm 
is  called  a  centralized  shortest  path  algorithm.   In  order  to 
apply  centralized  shortest  path  algorithms,  a  network  must 
establish  a  central  node  to  gather  information  concerning  the 
complete  network  topology  so  that  the  algorithm  can  be  executed 
at  the  central  node.   Many  authors  have  introduced  centralized 
shortest  path  algorithms  which  are  very  efficient  and  are  widely 
known  [4],  [5],  [9],  [10],  [14].    The  second  type  of  algorithm 
is  called  a  distributed  (or  decentralized)  shortest  path  algor- 
ithm.  Distributed  shortest  path  algorithms  do  not  require  a 
central  node  to  execute  the  algorithm;  instead,  they  are  executed 
locally  at  each  individual  node  using  only  the  limited  informa- 
tion received  from  its  neighboring  nodes.   As  compared  with 
centralized  shortest  path  algorithms,  distributed  shortest  path 
algorithms  are  less  efficient  and  less  well-known.   Authors  who 
have  introduced  distributed  shortest  path  algorithms  include: 
Abram  and  Rhodes  [1],  Boehm  and  Mobley  [3],  Fultz  and  Kleinrock 

[6],  McQuillan  [7],  Naylor  [8],  Segall,  Merlin  and  Gallager  [11] 
and  Tajibnapis  [12]. 


In  a  computer  network,  the  links  or  nodes  do  not  always 
function  properly.   They  often  change  from  a  state  of  function- 
ing to  a  state  of  breakdown.   Also,  the  length  of  a  link  may 
increase  or  decrease  from  time  to  time.   Therefore  it  is  necessary 
for  a  shortest  path  algorithm  to  adapt  to  such  topological  changes 
in  the  network.   In  our  discussion,  we  will  use  the  length  of 
links  to  indicate  whether  links  or  nodes  fail  or  recover. 
When  a  link  fails,  we  will  set  the  length  of  the  link  to  infinity; 
and  when  a  link  recovers,  we  will  change  the  link  length  from 
infinity  to  the  actual  length  of  the  link.   Similarly,  when  a 
node  fails,  we  will  set  the  length  of  all  links  associated  with 
the  node  to  infinity;  and,  when  a  node  recovers,  we  will  change 
the  associated  link  lengths  from  infinity  to  the  actual  length 
of  the  links.   Also,  the  term,  "adaptive  algorithm,"  is  used  to 
describe  an  algorithm  that  can  be  applied  to  a  network  where  the 
link  lengths  are  subject  to  periodical  changes. 

The  purpose  of  this  paper  is  to  present  an  adaptive 
distributed  algorithm  for  finding  shortest  path  from  a  source 
node  to  all  other  nodes  in  an  N-node  computer  network.   We  will 
present  this  algorithm,  first,  by  describing  a  basic  distributed 
algorithm  for  finding  shortest  paths  from  an  origin  to  all 
other  nodes  in  a  network  where  there  are  no  topological  changes 
(Algorithm  1) .   Then  we  will  present  the  modification  of 
Algorithm  1  that  is  applicable  to  topological  changes  in  networks 


(Algorithm  2) .  Also,  we  will  briefly  describe  how  these  algo- 
rithms can  be  applied  to  find  all  shortest  paths  between  every 
pair  of  nodes  in  the  network. 

The  advantages  of  Algorithm  2  are  as  follows: 

1)  The  algorithm  is  applicable  to  networks  with  directed  or 
undirected  links. 

2)  It  is  executed  locally  at  each  individual  node  using  only 
simple  information  obtained  from  the  neighborhing  nodes. 

3)  It  is  adaptive  to  topological  changes  such  as  changes  in 
link  lengths  and  the  failure  or  recovery  of  links  and 
nodes. 

4)  The  shortest  paths  obtained  at  the  termination  of  the 
algorithm  are  optimal. 

5)  The  algorithm  works  much  faster  than  other  available 

algorithms.   After  the  last  topological  change  is  made, 

3 
the  algorithm  requires  in  the  worst  case  a  total  of   N  /4 

3 
additions  and  N  /4   comparisons  to  obtain  all  optimal 

shortest  paths  from  a  source  node  to  all  other  nodes. 

However,  experience  has  shown  that  in  practice  the 

algorithm  terminates  much  faster  than  these  worst  bounds. 

6)  The  algorithm  is  easy  to  understand  and  easy  to  program. 


2.   ALGORITHM  1 

In  this  section  we  will  present  a  distributed  algorithm 
for  finding  the  shortest  paths  from  an  origin  to  all  other 
nodes  in  a  network  where  there  are  no  topological  changes.   In 
order  to  apply  the  algorithm  we  require  that: 

1)  Each  node  is  given  a  node  label:   node  1,  node  2,...,    node  N, 
where  node  1  is  the  origin. 

2)  Each  node  has  knowledge  of  a  set  of  neighboring  TO  nodes 
that  are  connected  to  node  I  by  directed  links  leading 
from  node  I  to  the  TO  nodes,  respectively.   And  the  set  of 
TO  nodes  are  further  separated  into  two  subsets:   HIGHER  TO 
nodes  and  LOWER  TO  nodes.   The  HIGHER  TO  nodes  are  the  TO 
nodes  whose  node  labels  are  larger  than  I,  and  the  LOWER  TO 
nodes  are  the  TO  nodes  whose  node  labels  are  smaller  than  I. 
For  example,  node  5  is  a  HIGHER  TO  node  of  node  2,  and  is 

a  LOWER  TO  node  of  node  9.   Each  node  I  also  has  knowledge 
of  the  lengths  of  the  directed  links  going  from  node  I  to 
the  TO  nodes. 

In  an  N-node  directed  network,  let 

I,  J,  L  =  1,2, ...,N,  be  the  nodes  of  the  network  where 

node  1  is  the  origin,  and  the  remain- 
ing nodes  are  numbered  arbitrarily, 
(I, J)   be  the  directed  link  connecting  node  I 
to  node  J, 
D(I,J)  >  0   be  the  length  of  link  (I, J), 


F(J:K),  J  =  2,3,... ,N   and   K  =  1,2,...   be  the  distance  of 
the  tentative  shortest  path  from  node  1  to  node  J 
at  the  current  cycle   K   of  the  algorithm, 
H(J:K)   be  the  node  that  precedes  node  J  on  a  tentative 

shortest  path  from  node  1  to  node  J  whose  distance 
is   F(J:K). 
More  notation  will  be  introduced  later  when  it  becomes  necessary. 

Algorithm  1  is  developed  from  the  following  principle. 
At  a  pre-arranged  time,  each  node  I,  I  =  2,3,...,N,  in  turn 
informs  each  of  its  TO  nodes  J  what  is  the  current  best  shortest 
distance  from  node  1  to  node  J  via  node  I.   This  procedure 
allows  the  shortest  distances  from  node  1  to  all  other  nodes  to 
be  updated  iteratively  in  order  to  reach  optimality.   The  follow- 
ing steps  can  be  taken  to  achieve  this  principle. 

1)  Divide  the  time  horizon  into   N   cycles  and  divide  each 
cycle  into   N   time  periods. 

2)  In  an  odd  cycle   K   at  time  period   T:  node  I  (where   I  =  T) 
sends  each  of  its  HIGHER  TO  nodes  J  the  current  best  shortest 
distance  from  node  1  to  node  J,  i.e.,  F(J:K). 

3)  In  an  even  cycle   K   at  time  period   T:  node  I  (where 

I  =  N-T+l)   sends  each  of  its  LOWER  TO  nodes  J  the  current 
best  shortest  distance  from  node  1  to  node  J,  i.e.,  F(J:K). 

4)  In  any  cycle  K,  any  node  I  is  not  to  send  its  TO  nodes 
the  current  shortest  distances   F(J:K)'s  unless 
F(J:K)  <  F(J:K-2) . 


A  proof  of  why  these  steps  will  determine  the  optimal  shortest 
distances  from  node  1  to  all  other  nodes  in  at  most  N-l  cycles 
can  be  found  in  [13],  [14,  pp.  52-71]. 

We  will  now  describe  the  distributed  algorithm  for 
finding  the  shortest  paths  from  node  1  to  all  other  nodes  in 
a  network  where  there  are  no  topological  changes.   Note  that 
the  following  Algorithm  1  is  to  be  executed  by  cycles  in  the 
order  K  =  0,1,2 ,... ,N-1. 

Algorithm  1 

Cycle  0    (to  initialize   F(J:-l)'s  and   F(J:0)*s) 

A.  Let   F(J:-1)  =  °°   for  all   J  =  1,2, ,N. 

B.  Let   F(1:0)  =  0,  F(J:0)  =  D(1,J)   and   H(J:0)  =  1   for  all 
J  =  2,3,. . . ,N. 

Cycle   K  =  1,3,...   (to  update   F(J:K)'s   for  pivot  node  I's 

HIGHER  TO  nodes  J) 

A.  Let   F(J:K)  =  F(J:K-1)   for   J  =  1,2,. . . ,N. 

B.  At  each  of  time   T,  T  =  1,2,..., N-l,  do  the  following: 

I.   (at  pivot  node  I,  where   I  =  T) 
Check  if   F(I:K)  <  F(I:K-2). 

a.  If   F(I:K)  _>  F(I:K-2),  do  nothing  and  go  to  time  T+l 

b.  If   F(I:K)  <  F(I:K-2),  compute   F(J:K)*   by 

F(J:K)*  =  F(I:K)  +  D(I,J)   for  each  of  node  I's 

HIGHER  TO  nodes  J;  and,  transmit  "F(J:K)*M  and  "I" 

to  the  corresponding  HIGHER  TO  nodes  J.   Go  to 

time  T+l. 
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II.   (at  each  of  HIGHER  TO  nodes  J  of  node  I) 

After  receiving  "F(J:K)*"  check  if   F(J:K)*  <  F(J:K). 

a.  If   F(J:K)*  >_   F(J:K),  do  nothing. 

b.  If   F(J:K)*  <  F(J:K),  F(J:K*)   becomes   F(J:K) 
and  I  becomes   H(J:K);  or,  F(J:K)*  <-  F(J:K) 
and   H(J:K)  <-    I. 

III.   (at  all  other  nodes) 
Do  nothing. 

Cycle   K  =  2,4,...  (to  update   F(J:K)'s  for  pivot  node  I's 

LOWER  TO  nodes  J) 

A.  Let   F(J:K)  =  F(J:K-1)   for   J  =  1,2,. ..,N. 

B.  At  each  of  time   T,  T  =  l,2,...,N-2,  do  the  following: 

I.   (at  pivot  node  I,  where   I  =  N-T+l) 
Check  if   F(I:K)  <  F(I:K-2). 

a.  If   F(I:K)  _>  F(I:K-2),  do  nothing  and  go  to 
next  time  T+l . 

b.  If   F(I:K)  <  F(I:K-2)   compute   F(J:K)*   by 
F(J:K)*  =  F(I:K)  +  D(I,J)   for  each  of  node  I's 
LOWER  TO  nodes  J,  and  transmit  "F(J:K)*M  and 
"I"  to  the  corresponding  LOWER  TO  node  J.   Go 
to  time  T+l. 


Start 


© 


© 


F(Jj-1)     -    -J    J    -    1 


,2, N~] 


F(1:0)    -    0,    F(J:0)    -    D(1,J)    and    H(J:0)     =    1    for    I    -    2,3,. ..,N 


K   =   0 


© 


Cycle   K   =   K+l 


F(J:K)    =    F(J:K-1),    J    =    1,2,.  .  .  ,N 


Yes 


~©T 


© 

No 


I    =    T 


F(I:K)     <    F(I:K-2)?     ~) 


XE 


Yes 


© 


J  I  {HIGHER  TO  nodes    of    1} 


) 


F(J:K) *    =    F(I:K}     +    D(I,J) 


© 

No 


Transmit    to   node   J 
*F(J:K) *"    and    "I" 


~C    F(J:K) «    <    F(J;K) ?       ) 
]  Yes   (b) 


F(J:K)    -    F(J:K)*,    H(J:K)    -     I 


No 


Xk    "    N-l?) 


r-i 


Yea 


Stop 


-{k    =    odd^- 


NO 


/  T  =    1,2,..., N-l    V*" 


© 


I    =   N-T+l 


©. 


No 


C    F(I::0     <    F(I:K-2)  ? 
Yes    (b) 


c 


J,    J      €  {LOWER  TO   nodes    of    1} 


F(J:K) *    =    F(I:K)     +    D(I,J) 


Transmit  to  node  J 
"F(J:K) *"  and  "I" 
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No 
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FIGURE  1.   A  flow  chart  of  Algorithm  1  for  finding  the 

shortest  paths  from  node  1  to  all  other  nodes 
The  circled  letters  are  the  step  numbers  of 
the  algorithms. 
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II.   (at  the  LOWER  TO  nodes  J  of  node  I) 

After  receiving   "F(J:K)*H,  check  if   F(J:K)*  <  F(J:K). 

a.  If   F  (J:K)  *  _>  F(J:K)  ,  do  nothing. 

b.  If   F(J:K)*  <  F(J:K),  let   F(J:K)  +    F(J:K)* 
and   H(J:K)  *■  I. 

III.   (at  all  other  nodes) 
Do  nothing. 

A  flow  chart  of  Algorithm  1  is  given  in  Figure  1. 

The  optimal  distances  of  the  shortest  paths  from  node  1 
to  all  other  nodes  are  determined  when   F(J:K)  =  F(J:K-1)   for 
all   J,  J  =  2,3,... ,N.   When  this  condition  occurs,  all  nodes  J 
will  cease  transmitting;  and,  this  condition  will  occur  no 
later  than  in  cycle  N-l.   We  will  use  an  example  to  illustrate 
how  Algorithm  1  determines  the  shortest  distances  from  node  1 
to  all  other  nodes  in  no  more  than  N-l  cycles.   Suppose  that 
in  a  10-node  network  the  shortest  path  from  node  1  to  node  9 
is  1-3-8-7-5-2-4-6-9.   Then  the  shortest  distance  from  node  1 
to  node  3  is  determined  in  cycle   K  =  1,  time  T  =  1   as 
F(3:l);  the  shortest  distance  to  node  8  is  determined  in 
K  =  1,  T  =  3   as   F(8:l) ;  to  node  7  in   K=2,T=3   as 
F(7:2);  to  node  5  in   K  =  2,  T  =  4   as   F(5:2);  to  node  2  in 
K  =  2,  T  =  6   as   F(2:2) ;  to  node  4  in   K  =  3,  T  =  2   as 
F(4:3);  to  node  6  in   K  =  3,  T  =  4   as   F(6:3);  and,  to  node  9 
in   K  =  3,  T  =  6   as   F(9:3).   As  a  matter  of  fact,  the  optimal 
distances  of  all  shortest  paths  from  node  1  to  other  nodes 
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that  include  a  sequence  of  successively-larger  nodes,  such  as 
1-3  and  1-3-8  in  our  example,  are  all  determined  in  cycle   K  =  1, 
In  cycle   K  =  2,  more  shortest  paths  are  determined  by  attaching 
a  sequence  of  successively-smaller  nodes,  such  as  7-5  and  7-5-2, 
to  the  nodes  of  the  shortest  paths  already  determined  in  cycle 
K  =  1.   In  cycle   K  =  3,  more  shortest  paths  are  determined  by 
attaching  a  sequence  of  successively-larger  nodes,  such  as 
4-6  and  4-6-9  in  our  example,  to  the  nodes  of  the  shortest  paths 
that  are  determined  in  cycle   K  =  2,...  .   Since  a  shortest 
path  can  include  at  most  N-l  successively-increasing  and 
successively-decreasing  sequences  of  nodes,  all  shortest  paths 
in  the  network  are  determined  in  no  more  than  N-l  cycles. 
Actually,  the  number  of  cycles  it  takes  the  algorithm  to 
determine  all  shortest  paths  is  equal  to  the  maximum  number  of 
successively- increasing  and  successively-decreasing  sequence 
of  nodes  in  any  of  the  shortest  paths  from  node  1  to  all  other 
nodes. 

When  the  algorithm  is  terminated,  say,  at  the  end  of 
cycle  N-l,  each  node  J,  J  =  2,3,...,N,  has  the  following 
solutions : 

1)  The  optimal  distance  of  the  shortest  path  from  node  1 
to  node  J;  i.e.,  F(J:N-1),  and 

2)  The  identity  of  the  node  that  proceeds  node  J  on  the 
shortest  path  from  node  1  to  node  J;  i.e.,  H(J:N-1). 
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However,  node  1  does  not  know  what  are  the  shortest  distances 
to  all  other  nodes  and  what  are  the  next  nodes  on  the  shortest 
paths  to  these  nodes.   If  desired,  node  1  can  obtain  this 
information  by  requiring  each  node  J  to  transmit   F(J:N-1) 
back  to  node  1.   If  link  (J, I)   exists  for  every  link  (I, J), 
node  J  can  easily  send   F(J:N-1)   back  to  node  1  via  the  inter- 
mediate nodes  on  the  shortest  path.   However,  if  link  (J, I)   does 
not  exist  for  every  link  (I.J),  node  J  has  to  find  a  new  path 
to  send  the   F(J:N-1)   back  to  node  1,  which  requires  sub- 
stantial effort. 

To  examine  the  efficiency  of  Algorithm  1,  recall  that 
there  are   N   nodes  in  the  network.   In  the  worst  case,  each 
of  cycles   K  =  1,3,...   requires  node  1  to  make   N-l   additions, 
comparisons  and  transmissions;  node  2  to  make   N-2   additions, 

comparisons  and  transmissions;  ...,  and  so  on.   The  sum  taken 

2 
over  nodes  is  approximately   N  /2   additions,  comparisons  and 

transmissions.   Similarly,  each  of  cycles   K  =  2,4,...   requires 

2 
a  total  of  approximately   N  /2   additions,  comparisons  and 

transmissions.   In  the  worst  case,  the  algorithm  has  to  be 

3 
executed  for   N-l   cycles,  therefore,  approximately   N  /2 

additions,  comparisons  and  transmissions  are  necessary  overall. 

However,  in  each  cycle  of  the  algorithm,  at  least  one  optimal 

shortest  distance  from  node  1  to  node  J  will  be  determined; 

and  the  determination  of  a  shortest  distance  will  save   N-l 

additions,  comparisons  and  transmissions  in  two  cycles  to  follow. 
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Consequently,  a  total  of   [ (N-l) 2/2]  [1  +  2  +  3 +  *  •  •  +  N-l]  =  N3/4 
additions,  comparisons  and  transmissions  are  saved  throughout 

the   N-l   cycles  of  the  algorithm.   Therefore,  in  the  worst 

3       3       3 
case,  Algorithm  1  requires  approximately   N  /2  -  N  /4  =  N  /4 

additions,  comparisons  and  transmissions.   In  the  best  case, 

Algorithm  1  can  determine  all  shortest  distances  in  one  cycle 

2 

using  only  N  /2   additions,  comparions  and  transmissions. 

Some  empirical  studies  have  found  that  algorithms  similar  to 
Algorithm  1  have  required  much  fewer  operations  than  the  worst 
computational  bound  of   N  /4   [14,  pp.  82-89]. 

The  relative  efficiency  of  Algorithm  1  appears  to  be 
quite  good  as  compared  with  other  available  algorithms.   For 
example,  the  worst  computational  bound  of  Algorithm  1  is  only 
one-half  of  the  worst  bounds  of  Abram  and  Rhodes1  [1]  and 
Segall,  Merlin  and  Gallager's  [11]  algorithms.   However,  in  most 
cases  Algorithm  1  will  work  much  faster  than  these  two  algorithms 
For  example,  in  a  10-node  network  where  the  shortest  paths  are: 

1-10-9-8-7-6-5-4-3-2,  Algorithm  1  will  require  two  cycles  using 

2 

approximately   N    steps;  whereas,  Abram  and  Rhodes'  algorithm 

2 
will  require  9  cycles  using  approximately   5N    steps.   When 

N   is  large  Algorithm  1  can  work  up  to  approximately   N/2 

times  faster  than  the  algorithm  of  Abram  and  Rhodes.   While 

Algorithm  1  is  expected  to  work  faster  than  Abram  and  Rhodes' 

algorithm  in  most  situations,  there  is  no  situation  that 

Algorithm  1  will  work  slower  than  their  algorithm. 
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Algorithm  1  represents  a  basic  algorithm  for  finding 
the  shortest  paths  from  node  1  to  all  other  nodes.   Algorithm  1 
can  be  modified  to  improve  its  use  and  efficiency.   For  example, 
node  1  can  be  excluded  from  execution  of  the  algorithm  after 
cycle  1. 

Algorithm  1  can  also  be  applied  to  find  the  shortest 
paths  from  any  node  to  all  other  nodes  in  the  network.   To  find 
the  shortest  paths  from  an  arbitrary  node  I  to  all  other  nodes, 
we  can  modify  the  algorithm  slightly  by  initializing   F(I:0)  =  0 
and   F(J:0)  =  <*>   for  all   I  ^  J   in  Step  B  of  Cycle  0,  and 
execute  cycles  1,2,...,  in  the  same  manner  as  Algorithm  1. 
Of  course  the   F(J:K)  *s  in  such  a  case  represent  the  shortest 
distances  from  node  I  to  node  J,  respectively.   However,  in  this 
case  the  algorithm  may  not  start  building  the  shortest  paths 
until  cycle   K  =  2,  time   T  =  N-T+l,  consequently,  the  algorithm 
should  be  executed  at  least  until  cycle   N,  time   T  =  N-i ,  to 
guarantee  an  optimal  solution.   For  simplicity,  when  the  origin 
is  not  node  1,  we  will  terminate  the  modified  Algorithm  1 
after  cycle  N  instead  of  N-l. 

The  above  modification  of  Algorithm  1  can  be  applied 
concurrently  to  find  the  shortest  paths  between  all  pairs  of 
nodes  in  N  cycles.   That  is,  in  each  cycle  K,  the  modified 
algorithm  can  be  executed  N   times  by  successively  taking 
node  L,  L  =  1,2,...,N,  as  the  source  node.   When  the  modified 
algorithm  is  applied  to  find  the  shortest  paths  between  all 
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4 
pairs  of  nodes  the  worst  computational  bound  is   N  /4   additions 

and  N  /4   comparisons.   However,  the  upper  bound  of  the 

3 
necessary  transmissions  between  the  nodes  is  only  N  /4,  because 

each  transmission  from  node  I  to  node  J  can  contain  the  shortest 
distances  from  all  N-l  source  nodes  to  node  J. 

Algorithm  1  can  also  be  modified  slightly  so  that  each 
pivot  node  I  does  not  have  to  transmit  messages  at  a  pre-arranged 
time.   This  can  be  done  as  follows:   In  an  odd  cycle,  node  I 
becomes  a  pivot  node  and  is  ready  to  transmit  whenever  it  has 
received  messages  from  all  of  its  neighboring  nodes  whose  node 
numbers  are  smaller  than  I;  and,  in  an  even  cycle,  node  I 
becomes  a  pivot  node  and  is  ready  to  transmit  whenever  it  has 
received  messages  from  all  of  its  neighboring  nodes  whose  node 
numbers  are  larger  than  I.   The  advantages  of  this  modification 
are  that  no  clock  is  necessary  at  each  node  and  that  the 
algorithm  can  be  executed  within  a  shorter  time  span.   This 
time  savings  occurs  because  instead  of  waiting  for  its  pre- 
arranged time  to  transmit  messages,  node  I  can  transmit  messages 
as  soon  as  it  becomes  a  pivot  node. 
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3.   ALGORITHM  2. 

In  this  section  we  will  present  Algorithm  2  for  finding 
the  shortest  paths  from  node  1  to  all  other  nodes  in  a  network 
which  is  adaptive  to  topological  changes.   Algorithm  2  is  a 
modification  of  Algorithm  1  and  uses  the  following  additional 
ideas : 

1)  When  the  link  length   D(I,J)   is  reduced,  the  algorithm 
updates   F(J:K)   using  the  reduced   D(I,J). 

2)  When  the  link  length   D(I,J)   is  increased  such  that  the 
value  of   F(J:K)   may  be  increased,  the  algorithm  discards 
all  tentative  shortest  paths  that  use  link  (I, J)  on  their 
paths,  and  then  rebuilds  new  paths. 

3)  When  the  last  link,  (I, J),  on  the  tentative  shortest  path 
to  node  I  fails,  the  algorithm  discards  all  tentative 
shortest  paths  that  use  link (I, J)  on  their  paths,  and  then 
rebuilds  new  paths. 

The  distributed  algorithm  for  finding  the  shortest 
paths  from  node  1  to  all  other  nodes  in  a  network  which  is 
adaptive  to  topological  changes  is  as  follows.   Note  that  in 
Algorithm  2,  unless  otherwise  specified,  D(I,J)   represents 
the  current  link  length  at  the  current  cycle  and  time  period, 
which  can  be  different  from  earlier  lengths  due  to  topological 
changes . 
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Algorithm  2 

Cycle  0  (to  initialize   F(J:0)'s) 

A.   Let   F(1:0)  =  0,  F(J:0)  =  D(lfJ)   and   H(J:0)  =  1   for 
J  =  2 , 3, . . . ,N. 

Cycle  K  =  1,3,...  (to  update   F(J:K)'s  for  the  pivot  node  I's 

HIGHER  TO  nodes  J) 

A.  Let   F(J:K)  =  F(J:K-1)   for  all   J  =  1,2,. ..,N. 

B.  At  each  of  time   T,  T  =  1,2,...,N,  do  the  following: 

I.   (at  pivot  node  I,  where   I  =  T) 

a.  Check  if   F(I:K)  =  °°. 

i.   If  yes,  go  to  Step  b. 
ii.   If  no,  go  to  Step  c. 

b.  Check  if   F(I:K)   was  set  to   °°   due  to  the 
message  "DISMANTLE"  received  since  the  last  time 
node  I  was  the  pivot  node  in  cycle  K-l,  time 

T  =  N-I+2. 

i.   If  no,  go  to  time  T+l. 
ii.   If  yes,  send  node  I's  HIGHER  TO  and  LOWER  TO 

node  messages:  "DISMANTLE"  and  "I".   Then 

go  to  time  T+l. 
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c.  Check  if  link  (L,I)  still  exists,  where  node  L, 
L  =  H(I:K) ,  is  the  node  from  which  node  I  has 
obtained  the  current   F(I:K). 

i.   If  yes,  go  to  Step  d. 
ii.   If  no,  let   F(I:K)  «-  °°   and  send  each  of 
node  I's  HIGHER  TO  and  LOWER  TO  nodes  J 
the  messages:   "DISMANTLE"  and  "I".   Then 
go  to  time  T+l. 

d.  For  each  of  node  I's  HIGHER  TO  nodes  J,  check 
if   D(I,J)   has  been  increased  since  cycle  K-2, 
time   T  =  1  +  1 . 

i.   If  no,  compute   F(J:K)*  =  F(I:K)  +  D(I,J) 

and  transmit  the  messages  "F(J:K)*"  and  "I" 

to  node  J. 

ii.   If  yes,  send  node  J  the  messages:  " DISMANTLE" 

and  "I".   Go  to  time  T+l  after  Step  d  is  done 

II.   (at  each  of  node  I's  TO  nodes  J) 

a.   If  node  J   receives  the  messages  "F(J:K)*"  and 

"I",  do  the  following: 

i.   Check  if   F(J:K)  <■    °°   due  to  the  message 

"DISMANTLE"  received  since  cycle  K-l, 

time  T  =  N-I+2 . 

If  yes,  do  nothing. 
If  no,  go  to  Step  ii. 
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ii.   Check  if   F(J:K)*  <  F(J:K). 
If  no,  do  nothing. 

If  yes,  let   F(J:K)  <-  F(J:K)*   and   H(J:K)  «-  I 
b.   If  node  J  receives  the  messages:  "DISMANTLE"  and 
"I",  check  if   I  =  H(J:K) . 

i.   If   I  ±   H(J:K),  do  nothing, 
ii.   If   I  =  H(J:K),  let   F(J:K)  *■  <*>,  and  do  not 
update   F(J:K)   until  after  node  J  becomes 
a  pivot  node. 
III.   (all  other  nodes) 
Do  nothing. 

Cycle  K  =  2,4....  (to  update   F(J:K)'s  for  the  pivot  node  I's 

LOWER  TO  nodes  J) 

A.  Let   F(J:K)  =  F(J:K-1)   for  all   J  =  1,2,. ..,N. 

B.  At  each  of  time  T,  T  =  1,2, ...,N,  do  the  following: 

I.   (at  pivot  node  I,  where   I  =  N-T+l) 

a.  Check  if   F(I:K)  =  °°. 

i.   If  yes,  go  to  Step  b. 
ii.   If  not,  go  to  Step  c. 

b.  Check  if   F(I:K)   was  set  to   °°   due  to  the 
message  "DISMANTLE"  received  since  the  last  time 
node  I  was  the  pivot  node  in  cycle  K-l,  time 

T  =  1+1. 

i.   If  no,  go  to  time   T+l . 

ii.   If  yes,  send  node  I's  LOWER  TO  and  HIGHER  TO 

nodes  messages:  "DISMANTLE"  and  "I".   Then  go 

to  time   T+l . 
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c.  Check  if  link  (L,I)  still  exists,  where  node  L, 

L  =  H(I:K),  is  the  node  from  which  node  I  has  obtained 

obtained  the  current   F(I:K). 
i.   If  yes,  go  to  Step  d. 

ii.   If  no,  let   F(I:K)  •*-  °°   and  send  each  of 
node  I's  LOWER  TO  and  HIGHER  TO  nodes  J 
the  messages:  "DISMANTLE"  and  "I".   Then 
go  to  time  T+l . 

d.  For  each  of  node  I's  LOWER  TO  nodes  J,  check 
if   D(I,J)   has  been  changed  since  cycle  K-2, 
time  T  =  1+1. 

i.   If   D(I,J)   has  not  change  or  if  it  has 

decreased,  compute   F(J:K)*  =  F(I:K)  +D(I,J) 
and  transmit  the  messages  "F(J:K)*"  and 
"I"  to  node  J. 
ii.   If   D(I,J)   has  increased,  send  node  J  the 
messages:  "DISMANTLE"  and  "I". 
Go  to  time  T+l  when  Step  d  is  done. 
II.   (at  each  of  node  I's  TO  nodes  J) 

a.   If  node  J  receives  the  messages  "F(J:K)*" 
and  "I",  do  the  following: 
i.   Check  if   F(J:K)  «-  °°   due  to  the  message 
"DISMANTLE"  received  since  cycle  K-l, 
time  T  =  1+1. 
If  yes,  do  nothing. 
If  no,  go  to  Step  ii. 
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ii.   Check  if   F(J:K)*  <  F(J:K). 
If  no,  do  nothing. 

If  yes,  let   F(J:K)  +■  F(J:K)*  and  H(J:K)  «-  I. 
b.   If  node  J  receives  the  messages:  "DISMANTLE" 
and  "I",  check  if   I  =  H(J:K). 

i.   If  I    ?   H(J:K),  do  nothing, 
ii.   If   I  =  H(J:K),  let   F(J:K)  «-  °°,  and  do 
not  update   F(J:K)   until  node  J  becomes 
a  pivot  node. 
III.   (at  all  other  nodes) 
Do  nothing. 
The  flow  chart  of  odd  cycles,  K  =  1,3,...,  of  Algorithm  2 
is  given  in  Figure  2.   Even  cycles  work  similarly. 

The  speed  with  which  Algorithm  2  will  terminate  depends 
on  when  the  last  topological  change  is  made.   In  the  worst  case, 
Algorithm  2  will  terminate  in   N   cycles  after  the  last  topo- 
logical change  is  made.   While  Algorithm  1  successively  builds 
shortest  paths,  Algorithm  2  simultaneously  builds  and  dismantles 
shortest  paths  whenever  such  operations  are  necessary.   When 
there  are  no  topological  changes,  Algorithm  2  works  similar  to 
Algorithm  1.   However,  when  there  are  topological  changes, 
Algorithm  2  will  adapt  to  the  changes  in  order  to  find  the 
shortest  paths  from  node  1  to  all  other  nodes.   We  will  examine 
how  Algorithm  2  adapts  to  topological  changes  in  cycle  K, 
time   T   as  follows.   Since  odd  cycles  work  similar  to  even 
cycles,  we  will  use  only  an  odd  cycle   K   in  our  discussion. 
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In  cycle   K   (K  =  1,3,...),  time  T  (T  =  1,2,...,N), 
a  pivot  node  I  (where   I  =  T)   is  associated  with  one  of  the 
following  conditions: 

1)  F(I:K)  =  ~. 

2)  F(I:K)  ±  °°  but  the  length  of  link  (I, J)  that  connects 
node  I  to  a  TO  node  J,  D(I,J)  has  not  changed  or  has 
reduced. 

3)  F(I:K)  ^   °°     but  the  link  length   D(I,J)   has  increased. 

4)  F(I:K)  f   °°      but  the  last  link  on  the  tentative  shortest 
path  to  node  I,  (L,I),  has  become  disconnected. 

We  will  examine  how  Algorithm  2  adapts  to  these  conditions 
as  follows: 

Case  1.   If   F(I:K)  =  °°,  there  is  no  use  trying  to  update 

F(J:K)   by   F(J:K)*  =  F(I:K)  +  D(I,J),  where  node  J 
is  a  HIGHER  TO  node  of  node  I.   Steps  B.I.a.i  and 
I.B.b  prevent   F(J:K) 's  from  being  updated  when 
F(I:K)  =  °o.   However,  F(J:K)  =  °°   may  also  indicate 
that  a  tentative  shortest  path  to  node  J  existed 
previously,  but  has  to  be  dismantled  because  the 
last  link  of  the  path,  (L,I),  either  has  failed  or 
has  increased  its  length  since  last  time  node  I 
was  a  pivot  node.   In  such  a  case,  F(I:K)  =  °° 
should  not  be  used  to  update   F(J:K) 's,  and  node  I 
must  notify  nodes  J  of  this  fact  so  that  they  can 
dismantle  their  tentative  shortest  paths  in  case 
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these  paths  go  through  node  I  and  use  link  (L,I) 
on  their  paths.   Steps  B.I.b.ii  of  the  algorithm 
serves  to  warn  its  TO  nodes  J  of  the  fact  that  node  I 
has  just  dismantled  its  tentative  shortest  path  and 
that  nodes  J  should  also  dismantle  their  tentative 
shortest  paths  in  case  these  paths  go  through  node  I 
and  use  link  (L,I)  on  their  paths. 

Case  2.   If   F(I:K)  ^  °°   and  the  link  length   D(I,J)   has 
been  reduced,  F(J:K)   must  be  updated  by 
F(J:K)*  =  F(I:K)  +  D(I,J)   using  the  smaller  D(I,J). 
The  same  process  must  also  be  performed  to  update 
F(J:K)   when  there  is  no  change  in   D(I,J)   but 
F(I:K)  <  F(I:K-2) .   In  order  to  simplify  the  descrip- 
tion of  the  algorithm,  Step  B.I.d.i  of  the  algorthm 
updates   F(J:K)   whenever   F(I:K)  ^  °°   and   D(I,J) 
is  the  same  or  has  become  smaller. 

Case  3.   When  the  length  of  link  (I, J)  on  a  tentative  shortest 
path  is  increased,  the  current  tentative  shortest 
path  to  node  J  may  have  to  be  replaced.   This  is 
necessary  because  the  increase  in  link  length 
D(I,J)   may  cause  the  current  shortest  path  to  become 
longer  than  other  available  paths.   Step  B.I.d  of 
the  algorithm  first  examines  if   D(I,J)   has  been 
increased  since  cycle  K-2,  time  T  =  1+1;  and,  if  so, 
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Step  B.I.d.ii  sends  node  J  the  messages:  "DISMANTLE" 
and  "I".   When  node  J  receives  these  messages,  Step 
B.II.b  of  the  algorithm  examines  if  the  tentative 
shortest  path  to  node  J  goes  through  node  I  and 
uses  link  (I, J),  by  checking  to  see  if   I  =  H(J:K). 
If   I  =  H(J:K),  the  tentative  shortest  path  to  node  J 
goes  through  node  I  and  uses  link  (I,J)  on  its  path. 
In  such  a  case  node  J  dismantles  the  tentative 
shortest  path  by  setting   F(J:K)  =  °°   in  Step  B.ll.b.ii, 
and  informs  its  own  TO  nodes  to  do  the  same  the  next 
time  node  J  becomes  a  pivot  node.   In  the  meantime, 
in  order  to  help  node  J  remember  the  fact  that   D(I,J) 
has  increased,  Step  B.I.b.ii  and  B.II.a.i  of  the 
algorithm  keep   F(J:K)  =  °°   from  being  updated  until 
node  J  has  its  opportunity  to  send  its  own  TO  nodes 
the  messages  "DISMANTLE"  and  "J"  when   J   becomes  a 
pivot  node.   Thus,  with  the  help  of  Step  B.I.b.ii, 
the  dismantling  process  goes  on  in  the  rest  of  cycle 
and  in  later  cycles,  until  all  tentative  shortest 
paths  that  use  link (L, I)  are  dismantled. 

Case  4.   When  link  (L,I) ,  the  last  link  on  the  tentative  shortest 
path  to  node  J,  becomes  disconnected,  node  L  cannot 
communicate  directly  with  node  J  any  longer.   There- 
fore, node  J  must  find  out  about  the  disruption  of 
link  (L,I)  by  itself.   Prior  to  updating   F(J:K)*s 
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for  node  I's  HIGHER  TO  nodes  J,  Step  B.I.c  of  the 
algorithm  checks  to  see  if  the  link (L, I)  still 
exists.  If  link  (L,I)  still  exists,  the  algorithm 
will  proceed  to  carry  out  other  necessary  steps  of 
the  algorithm.   However,  if  link  (L,I)  no  longer 
exists,  Step  B.I.c. ii  of  the  algorithm  will  set 
F(I:K)  =  °°   and  send  node  I's  TO  nodes  the  message. 
"DISMANTLE"  and  "I".   Then,  if  necessary,  nodes  J 
will  dismantle  the  tentative  shortest  paths  to 
nodes  J  by  setting   F(J:K)  =  °°,  and  send  the  messages 
further  on. 

A  link  in  the  network  can  be  associated  only  with 
on  of  three  topological  changes:   a  decrease  in  link  length 
(including  recovery  of  the  link) ,  an  increase  in  link  length, 
and  breakdown  of  the  link.   As  described  above,  Algorithm  2 
can  adapt  to  any  such  changes  at  any  links,  therefore,  the 
algorithm  is  adaptive  to  networks  where  there  are  topological 
changes. 

The  speed  with  which  Algorithm  2  will  terminate  depends 
on  when  the  last  topological  change  is  made  in  the  network. 
When  there  is  no  topological  change  in  the  network,  Algorithm  2 
has  an  efficiency  similar  to  that  of  Algorithm  1.   That  is, 

in  the  worst  case,  Algorithm  2  will  terminated  in  N-l  cycles 

3 

using  approximately   N  /4   steps  and,  in  the  best  case,  will 

2 

terminate  in  two  cycles  using   N    steps. 
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When  there  are  topological  changes  in  cycle   K   of 
the  algorithm,  the  algorithm  will  require  in  the  best  case 
practically  no  additional  effort  in  order  to  find  all  shortest 
paths,  simply  because  these  topological  changes  have  no  effect 
on  the  current  solutions.   In  the  worst  case,  when  there  are 
link  recoveries  and  link  length  reductions  in  cycle   K,  all 
shortest  paths  have  to  be  determined  anew  so  that  the  algorithm 
has  to  be  executed  N-l   more  cycles.   Also  in  the  worst  case, 
when  there  are  link  failures  or  increases  in  link  lengths  in 
cycle  K,  the  algorithm  will  require  up  to   N-l   additional 
cycles  to  dismantle  all  relevant  shortest  paths  and  up  to   N-l 
additional  cycles  to  update  the  new  shortest  paths.   However, 
the  update  cycles  start  only  about  one  cycle  behind  the  dis- 
mantling cycles,  and  then  are  executed  simultaneously  with  the 
dismantling  cycles.   Therefore,  when  there  are  link  failures 
or  increases  in  link  lengths  in  the  network,  the  algorithm 
requires,  in  the  worst  case,  an  additional   N   cycles  to 
determine  all  new  shortest  paths.   Overall  speaking,  after  the 
last  topological  change  is  made,  Algorithm  2  determines  all 

shortest  paths  from  node  1  to  all  other  nodes  in   N   cycles 

3 
using  up  to   N  /4   steps.   While  the  computational  upper  bound 

of  Algorithm  2  is  one-half  or  less  of  Abram  and  Rhodes'  and 

Segall,  Merlin  and  Gallager's  algorithms,  as  illustrated  in 

example  in  the  last  part  of  Section  2,  there  are  situations 

where  Algorithm  2  can  run  up  to   N/2   times  faster  than  these 

two  algorithms. 
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Algorithm  2  represents  a  basic  adaptive  algorithm 
for  finding  the  shortest  path  from  node  1  to  all  other  nodes. 
Algorithm  2  can  be  modified  to  improve  its  use  and  efficiency. 
For  example,  each  node  I  can  store  in  its  memory  all   F(I:K)'s 
that  node  I  has  received  from  all  other  nodes.   Then  when 
node  I  has  to  dismantle  the  best  tentative  shortest  path  to 
node  I,  it  can  replace  this  path  immediately  with  the  path 
that  has  the  next  smallest   F(I:K).   Of  course,  node  I  can 
compute   F(J:K)'s  and  send  them  to  its  TO  nodes  J  immediately 
after  the  best   F(J:K)   is  replaced  by  the  second-best   F(J:K). 
However,  the  messages  "DISMANTLE"  and  "I"  still  should  be  sent 
to  nodes  J  so  that,  if  necessary,  they  can  dismantle  their 
best  tentative  shortest  paths  and  replace  them  with  the  second- 
best  tentative  shortest  paths. 

Algorithm  2  can  also  be  modified  so  that  it  can  be 
used  to  find  the  shortest  paths  between  all  pairs  of  nodes 
in  N  cycles.   Or,  it  can  be  modified  so  that  pivot  nodes  I 
do  not  have  to  transmit  messages  at  a  prearranged  time.   We 
have  already  discussed  how  to  make  these  two  modifications 
with  respect  to  Algorithm  1. 
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